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ABSTRACT 

We present long-slit spectroscopy, continuum and [OIIIJA5007 imaging data ob- 
tained with the Very Large Telescope and the Gran Telescopio Canarias of the type 2 
quasar SDSS JO 123+00 at z =0.399. The quasar lies in a complex, gas-rich environ- 
ment. It appears to be physically connected by a tidal bridge to another galaxy at a 
projected distance of ~100 kpc, which suggests this is an interacting system. Ionized 
gas is detected to a distance of at least ~133 kpc from the nucleus. The nebula has a 
total extension of ~180 kpc. This is one of the largest ionized nebulae ever detected 
associated with an active galaxy. Based on the environmental properties, we propose 
that the origin of the nebula is tidal debris from a galactic encounter, which could 
as well be the triggering mechanism of the nuclear activity. SDSS J0123+00 demon- 
strates that giant, luminous ionized nebulae can exist associated with type 2 quasars 
of low radio luminosities, contrary to expectations based on type 1 quasar studies. 

Key words: (galaxies:) quasarsiemission lines; (galaxies:) quasars individual: SDSS 
J012341.47+004435.9; galaxies:ISM; 



1 INTRODUCTION 

Previous works suggest that the most luminous extended 
emission line regions (EELRs) around type 1 quasars 
(QSOs) at redshift z <0.5 exist preferentially associ- 
ated with steep-spectrum radio-loud QSOs with luminous 
[OIII]A5007 nuclear emission and low broad line region 
metallicities {Z <O.6Z0). These EELRs have low metallic- 
ities as well and seem especially prevalent in QSOs show- 
ing signs of strong interaction. The origin of these nebulae 
is controversial, whether cold accretion of intergalactic gas, 
tidal debris from a galactic merger or remnants of galactic 
superwinds (for a review see Fu & Stockton 2009; Stockton, 
Fu & Canalizo 2006). 

Very little is known about the existence of EELRs as- 
sociated with type 2 quasars (e.g. Humphrey et al. 2009; 

* Based on observations carried out at the European Southern 
Observatory (Paranal, Chile) with FORS2 on VLT-UTl (pro- 
gramme 69.A-0123(A)) and at the Observatory Roque de los 
Muchachos (La Palma, Spain) with OSIRIS on GTC (programme 
GTC27-09A). 



Gandhi, Fabian & Crawford 1996). These objects are unique 
laboratories to investigate the existence and properties of 
quasar EELR. The fortuitous occupation of the active galac- 
tic nucleus (AGN) acts like a "natural coronograph" , allow- 
ing a detailed study of many properties of the surrounding 
medium, without the problems associated with the bright 
quasar point spread function. On the other hand, due to 
their much lower radio luminosity, distortions imprinted by 
the radio activity on the properties of the underlying nebu- 
lae are less important than in radio-loud objects. 

For these reasons, we are undertaking a project based 
on imaging, long-slit and integral field data obtained with 
F0RS2/VLT, the integral field spectrograph PMAS on the 
3.5m telescope at Calar Alto Observatory (Humphrey et al. 
2009), and the tunable filter OSIRIS on the GTC. These are 
being used to investigate the existence of EELR associated 
with Sloan Digital Sky Survey (SDSS) type 2 quasars at 
z ^0.3-0.4 and characterize their morphological, ionization 
and kinematic properties. 

In this paper we present results on the type 2 quasar 
SDSS J012341.47+004435.9 (SDSS J0123+00 hereafter) at 
z =0.399 (Zakamska et al. 2003). The object was detected 
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by the NVSS survey with a flux density of 11.7±0.6 mJy 
(Condon et al. 1998). The radio emission is spatially unre- 
solved. Assuming a radio spectral index of a = 0.7 (where 

oc z^~"), the expected radio luminosity density at 5 
GHz is 5x10^° erg s~"^ Hz~-^ sr~-^. This value suggests that 
the quasar is "radio-quiet" or at most "radio-intermediate" 
(Miller, Peacock & Mead 1990). 

We assume Qa=0.7, Qm=0.3, Ho=71 km s~^ Mpc~^ 
In this Cosmology, 1" corresponds to 5.33 kpc at ;2=0.399. 



2 OBSERVATIONS AND ANALYSIS 

The VLT observations were carried out on the 8th and 9th 
Sept 2008 using F0RS2 on UTl. The 691_55+69 filter was 
used to obtain the intermediate band image on the 1st night. 
It has a spectral window 6635-7185 A which covers both 
continuum and redshifted H/3 and [GUI] AA4959, 5007. The 
pixel scale is 0.25 arcsec pix~^. A diamond pattern mosaic 
of 4 images with individual exposures of 200 seconds was 
obtained, shifting each image by 10^' in RA or Dec relative 
to the previous one. This was done to minimize the im- 
pact of bad pixels and ghosts on the CCD. The images were 
combined, with simultaneous cosmic ray removal after bias 
subtraction. The seeing size measured from the final image 
is FWHM=0.95±0.05'' for the 1st night. 

Long slit spectroscopy was obtained at position an- 
gles (PA) PA60 (first night) and PA 76 (second night), 
aligned with extended diffuse emission or nearby poten- 
tially companion objects revealed by the images (Fig. 1). 
The 600RI+19 grism and the GG435+81 order sorting fil- 
ter were used with a 1^^ wide slit. The useful spectral range 
was ^^5030-8250 A corresponding to rest frame 3595-5901 
A. The object was moved along the slit between two differ- 
ent positions for a more accurate sky subtraction. The data 
were reduced following standard procedures with STAR- 
LINK and IRAF packages. The spectra were debiased, flat 
fielded, wavelength (A) calibrated, background subtracted, 
combined with simultaneous cosmic ray removal, flux cali- 
brated and corrected for Galactic extinction. Geometric dis- 
tortion was found to be negligible. 

The nights were photometric. Comparison of several 
spectrophotometric standard stars taken with a b" slit dur- 
ing the run gave a flux calibration accuracy of 5% over the 
entire spectral range. To perform the spectroscopic analysis, 
the line profiles were fitted with Gaussian functions. The 
FWHM values were corrected for instrumental broadening 
(IP=7.2zb0.2 A as measured from the sky emission lines). 
Point sources did not fill the slit during the second night 
(seeing FWHM=0.65ili0.05^^) and this introduces uncertain- 
ties on the kinematic analysis. For this reason, the PA60 
spectrum was used to characterize the kinematic properties 
at the different spatial regions identified in the spectra. 

GTC images were obtained with the OSIRIS camera 
on 17th Oct 2009 using the order sorting filter 694/44. The 
2x2 detector binning provided a spatial scale of 0.254 arcsec 
pix~^. A continuum image and a continuum-|-[OIII]A5007 
image were produced. Each one consisted of three 200 sec- 
ond exposures obtained through the tunable filters with 
FWHM=16 A and central A (Ac) at the optical axis of 7019 
A and 7069 A respectively. Due to the nature of the instru- 
ment Ac changes with radial distance from the optical axis. 



Ac values were 7002.9 and 7052.7 A for the on-band and 
off-band images respectively at the location of the target, 
corresponding to rest-frame [OIII]A5007 and the adjacent 
continuum. The average seeing was 0.8=b0.1''. The observa- 
tions were reduced and flux calibrated with the standard 
star G158-100 using IRAF routines. Tunable filter profiles 
were generated at Ac corresponding to the spatial location 
of G158-100 on the CCD and applied to its absolute flux 
spectrum to obtain the flux calibration factors. The indi- 
vidual target exposures for each filter were sky subtracted, 
spatially aligned and combined. 



3 RESULTS 

3.1 Imaging data 

The VLT intermediate band image is shown in Fig. 1 (top 
panel). The bottom panel shows the GTC pure [GUI] 5007 
(green) and continuum (red) composite image. The origi- 
nal [OIII]+continuum image was continuum subtracted, and 
therefore the green colour shows the pure gas emission. The 
central elliptical galaxy presents distorted contours due to a 
detached "blob" which is more evident in the GTC images 
and emits both continuum and strong line emission. Low 
surface brightness diffuse features and knots are detected 
well outside the optical size of the galaxy. A long filament 
extends roughly in the E-W direction, to a maximum pro- 
jected distance from the nuclear region of r r^lA" (75 kpc) 
and a total extension of 25^^ (133 kpc) (throughout the text 
distances and sizes are projected values). 

The filament is clearly detected in [OIII] extending to- 
wards and overlapping with galaxy "Gl", which emits not 
only continuum but also line emission since it is clearly de- 
tected in the pure [OIII] GTC image. Therefore it is at sim- 
ilar z as the quasar. Both objects seem to be physically 
connected by the filament mentioned above, which is rem- 
iniscent of a tidal bridge. "G2", which is detected both in 
continuum and [OIII], also lies at a similar z and is located 
roughly along the direction of the filament. Other galaxies 
are found near the quasar in projection, but we cannot con- 
firm that any of them are at similar z. 

The VLT and GTC images show that SDSS J0123+00 
lies in a rich gaseous environment and suggest that it is a 
member of an interacting system. 

3.2 Long slit spectra. 

The 2-dim spectra at PAGO and PA76 including the 
[OIII] AA4959, 5007 doublet are presented in Fig. 2. The low 
surface brightness features revealed by the images (Fig. 1) 
along the slit direction emit a rich emission line spectrum. 
Continuum emission is not detected from these features. 
The spectra reveal a total extension of the ionized nebula of 
12.5'' (67 kpc) and -34'' (180 kpc) along PAGO and PA76 
respectively. The maximum radial extent from the nucleus 
is r —25" (133 kpc) measured along PA 76 towards the East. 
The filament revealed by the images is part of this structure. 

The visual inspection of the 2-dim spectra reveal 3 dis- 
tinct spatial regions where the gas is characterized by differ- 
ent kinematic and ionization properties, which are recognis- 
able at both PA. We have named them the central region. 
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Figure 1. Top: FORS2-VLT intermediate band image of SDSS 
J0123+00 containing continuum, R(3 and [OIII] AA4959,5007 emis- 
sion. The field size is 48^^x31.5^^ or 254 kpc x 168 kpc. The con- 
tours (arbitrary units) highhght the low surface brightness fea- 
tures. The slit positions used for the VLT spectra are indicated. 
Bottom: GTC continuum subtracted [OIII] 5007 image (green) 
and continuum image (red), color coded in the region of interest. 
The images demonstrate that "Gl" and "G2" are at the same z 
as the QSO. The gaseous filament is reminiscent of a tidal bridge 
connecting the "Gl" galaxy and the quasar. The images suggest 
that SDSS J0123+00 is a member of an interacting system. 



the intermediate region and the low surface brightness neb- 
ula (LSBN hereafter). We have analysed them separately. 
Since it is possible that within each region there is a spa- 
tial variation of the kinematic and ionization properties, we 
have defined within each region one or more apertures (Ap. 
hereafter) as specified below. 

In order to characterize the ionization properties of the 
gas at each spatial location, we plot the line ratios for the dif- 
ferent apertures in the [0III]/H/3 vs. [Oil] A3727/[OIII] A5007 
diagnostic diagram (Fig. 3), which involves the strongest 
emission lines detected in most apertures. Only apertures for 
which at least two emission lines are detected are shown. For 
comparison, the standard photoionization model sequence 
often applied to low redshift type 2 active galaxies is also 
shown (see Villar-Martm et al. 2008 for a discussion of this 
model sequence applied to type 2 quasars). The models as- 
sume solar metallicity, gas density n=100 cm~^ and power- 
law index a=-1.5. The ionization parameter U varies along 
the sequence. A subsample of HII galaxies extracted from 
the catalogue by Terlevich et al. (1991) is also shown. 

• r < 1" . The central region (Ap. 4 in Fig. 2 and 3). 
This is the region in which the continuum emission is 
brightest. The continuum spatial profile is strongly peaked 
and is dominated by an unresolved central component 
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Figure 2. 2-dim FORS2-VLT spectra of the [OIII] doublet at 
PA60 (top) and PA76 (bottom). The spectra are shown in two 
different color scales to highlight the low surface brightness re- 
gions (left) or the regions closer to the nucleus (right). The aper- 
tures used for the spectroscopic analysis (see text) are shown with 
numbers. Notice the striking extension ~180 kpc of the knotty, 
extended nebula along PA76. The extension to the East encom- 
passes the filament in Fig. 1. 



(FWHM=0.92±0.03'', PAGO), although diffuse continuum 
emission is also detected in excess of the PSF wings up 
to a total extension of ^^A" or ^21 kpc. The central re- 
gion emits a rich emission line spectrum, typical of AGN 
(Fig. 3 dark blue solid circle), including strong HeIIA4686. 
The FWHM([OIII]) value measured from the 1-dim spec- 
trum is 430±30 km s"\ 

We have fitted the continuum in Ap. 4 stellar population 
models. The nebular continuum was first removed. Its con- 
tribution, 6% of the total continuum flux at 4020 A {z =0), 
was estimated from the H/3 flux. A reasonable range of pa- 
rameters (gas metallicities, temperature and density) and 
no reddening (as implied by the H/3/H7=0.48±0.07 value) 
were considered. The fit (Fig. 4) to the nebular continuum 
subtracted spectrum implies the presence of an old stellar 
population and an excess of blue light which cannot be ac- 
counted for with only old stars. If it is due to young stars, 
the best fit is obtained with an old stellar population of >5- 
10 Gyr of age that contributes with ^80% of the continuum 
at 4020 A {z =0) plus a young stellar population (age^ 
10 Myr) which contributes ^-^20%. However, based on these 
data alone we cannot distinguish if this blue excess is alter- 
natively due to scattered nuclear radiation (the polarization 
level measured by Zakamska et al. (2006) was (2.7=1=2.2%). 

• 1 " < r < 2.5" . The intermediate region (Ap. 3 and 
5 Fig. 2 and 3). It extends at both sides of the central re- 
gion, up to a distance of 2.0^^ to the E and r^2.Q" to the 
W respectively from the continuum centroid (PAGO). The 
location of Ap. 3 and 5 in the diagnostic diagram (Fig. 3, 
cyan solid square and triangle) shows that they are predom- 
inantly photoionized by the AGN. This is also confirmed by 
the presence of strong Hell emission. 

An interesting feature within the intermediate region. 
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Figure 3. Diagnostic diagram and the gas ionization properties. 
Different symbols are used for different apertures. PA60: Solid 
symbols. PA76: Hollow symbols. Blue and cyan symbols corre- 
spond to the central and intermediate regions apertures, while 
red symbols correspond to apertures across the LSBN. The small 
magenta symbols are HII galaxies from Terlevich et al. (1992) 
catalogue. The solid black line represents standard AGN pho- 
toionization models (see text). Numbers -1 and -2 indicate the 
log(U) values corresponding to the nearest models. 




Figure 4. Nuclear spectrum (black) (the nebular continuum has 
been removed) and the best fit stellar population model (red) 
obtained with the code STARLIGHT (Gid Fernandes et al. 2005). 
An excess of blue light above the old stellar population (yellow) 
is implied by the analysis. We do not know whether this excess is 
due to young stars (blue) and/or scattered AGN light. 



which is clearly seen in the GTC images (Fig. 1) is 
a detached blob, which looks rather compact (Ap. 5, 
Fig. 2). FWHM values along PAGO (1.15±0.05'') and PA76 
(l.lzbO.l^^) implying a physical size of a few kpc. Its centroid 
is located at ^1.3^^ from the continuum centroid. Faint con- 
tinuum emission is detected. The lines are blueshifted rel- 
ative to the rest of the ionized gas (Fig. 2) and show very 
asymmetric line profiles which reveal complex kinematics. 

• r >y^ The low surface brightness nebula (LSBN) has 
a total extension of ^180 kpc, as mentioned above (Fig. 2), 
and overlaps towards the East with the filament or tidal 
bridge identified in the images (Fig. 1). The fact that the 
FIRST and NVSS fluxes are consistent with each other 
shows that no significant radio emission at 1.4 GHz is present 
on scales >5^^ or 27 kpc and therefore any jet present is prob- 
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Figure 5. LSBN spectrum integrated from Ap. 6 to 9 along PA76. 
Flux in units of lO"-*^^ erg s"-*^ cm~^ A"-*^. 



ably much smaller than the nebula. The spectrum integrated 
across Ap. 6 to 9 along PA 76 is shown in Fig. 5. 

A study of the emission line kinematics of the LSBN along 
PAGO shows that they are quiescent: it emits very narrow, 
unresolved emission lines with FWHM<100 to <140 km s~^ 
depending on spatial location, and the rest-frame velocity 
shifts relative to the nucleus are also small (AV ^200 km 
s~^). The lines are narrower than the IP measured from the 
sky lines at most spatial locations along PA 76, confirming 
the general narrowness of the lines also in this direction. 

In order to investigate the spatial variation of the ioniza- 
tion properties of the LSBN, several apertures where consid- 
ered along the slit. Since the LSBN is knotty and inhomoge- 
neous, we selected apertures that isolate individual features 
when possible (Fig. 2). We have plotted the location of the 
apertures along PA60 and PA76 in the diagnostic diagram 
(Fig. 3). It is important to note that the standard AGN 
photoionization models, although successful at reproducing 
many spectral properties of type 2 active galaxies in gen- 
eral and type 2 quasars in particular, cannot explain the 
large scatter in the diagnostic diagrams presented by this 
class of objects (e.g. Robinson et al. 1987), with a substan- 
tial fraction overlapping with HII galaxies. The locations of 
the LSBN apertures of SDSS J0123+00 in the diagrams are 
consistent with the scatter presented by the general pop- 
ulation of type 2 active galaxies. Although a range of gas 
and/or AGN continuum properties has usually been claimed 
as responsible for the scatter, Villar-Martm et al. (2008) 
proposed that an alternative possibility is that young stars 
contribute to the gas excitation, in addition to the AGN. 
We explore this possibility further for SDSS J0123+00. 

Fig. 3 shows that Ap. 2 and 6 line ratios along PA60, 
which are closer to the AGN, are consistent with pure AGN 
photoionization, while moving outwards (Ap. 1, at r ^6.6^^ 
or 35 kpc), the line ratios are consistent with HII star form- 
ing galaxies. A similar trend is found along PA 76: Ap. 6 
presents an intermediate spectrum, with line ratios between 
the standard AGN models and HII galaxies (Hell is also de- 
tected). All other apertures (r >4.5^^ or 24 kpc, most clearly 
Ap. 2, 7, 8, 9) overlap with the HII galaxy zone. 

Therefore, although with the available information it is 
not possible to discard that the entire nebula is photoionized 
by the AGN, we find that the spatial variation of the LSBN 
line ratios can be explained if the inner regions are photoion- 
ized by the AGN, while moving outwards stellar photoion- 
ization becomes dominant. Line ratios based on [OI]A6300, 
Ha, [NII]]A6583, [SII]AA6716,6732 would help to discrimi- 
nate between the two scenarios. 
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In the stellar photoionization scenario, we can estimate 
the gas oxygen abundance in regions purely ionized by stars, 
using the R23 method (Pilyugin k, Thuan 2005). The spec- 
trum integrated across the LSBN has H7/H/3=0.43=b0.07, 
consistent within the errors with the theoretical case B re- 
combination value, and thus we assume no line reddening. To 
make sure that AGN photoionization is relatively negligible, 
we have used the spectra of those apertures along PA 76 (e.g. 
Ap. 7,8,9) which are unambiguously placed among the HII 
galaxies and far from the standard AGN models in Fig. 3, 
taking the errorbars into account. 

The individual spectra, as well as the integrated spec- 
tra across all these apertures, are too noisy to set useful 
constraints on the [OIII]A4363 flux and Tg. This introduces 
uncertainties on the abundance determination, due to the 
bimodal shape of the 12+log(0/H) vs. log(R23) function 
(Pilyugin & Thuan 2005). The i?23 log values measured for 
Ap. 7, 8 an 9 along PA76 are 1.01±0.06, 0.96±0.09 and 
1.10±0.06 respectively and imply abundances well below so- 
lar, in the range 8% to 42% solar taking errors and uncer- 
tainties into account, and depending on the aperture. Sim- 
ilarly abundances are often found in compact HII galaxies 
(e.g. Guzman et al. 1997, Terlevich et al. 1991). For com- 
parison, we have used the SDSS spectrum to constrain the 
abundance within the inner AGN photoionized region using 
the method discussed in Humphrey et al. (2008). The re- 
sulting oxygen abundance must be >25% solar in order to 
reproduce the optical line ratios. 

The narrowness of the spectral lines compared with the IP 
(see §2), implies that the spatial extent in the direction per- 
pendicular to the slit is <l" or <5.3 kpc. For the brightest 
knots isolated along the slit, FWHM^l.lzbO.l'' correspond- 
ing to 4.7±0.6 kpc, taking the seeing into account. Sizes 
<5 kpc are also measured for the UV rest frame continuum 
knots identified in the HST images of Zakamska et al. (2006). 

The H/3 luminosity for Ap. 7,8,9 is in the range (0.4- 
1.0) X 10^° erg s-^ If the gas is photoionized by stars, star 
forming rates (SFR) --0.09-0.23 M© yr"^ are implied (Ken- 
nicutt 1998; given the large EW of the lines, stellar absorp- 
tion of H/3 is expected to be negligible). Considering the 
small knot sizes, slit flux losses could be <30%. This un- 
certainty does not affect our conclusions. Similar SFRs and 
sizes are found for some nearby compact HII galaxies (Telles, 
Munoz-Tunon, Tenorio-Tagle 2001, Guzman et al. 1997). 



4 DISCUSSION AND CONCLUSIONS 

SDSS0123+00 at z =0.399 is associated with a giant, lumi- 
nous ionized nebula, which extends for >180 kpc. This is 
one of the largest EELRs ever detected associated with an 
AGN. The EELR L^oiii] -^3x10^^ erg s"^ is similar to that 
of very luminous EELRs associated with type 1 quasars. 
Therefore, this object demonstrates that very luminous and 
giant EELR can also be associated with low radio lumi- 
nosity quasars, contrary to expectations based on type 1 
quasar studies (see §1). The low EELR abundances inferred 
(<O.5Z0) and the clear signs of interactions are features of- 
ten found in type 1 quasars with luminous EELR. 

SDSS0123+00 lies in a gas rich environment and it 
seems to be a member of an interacting system. The quasar 
appears to be physically connected by a tidal bridge to a 



companion galaxy at '^lOO kpc in projection. The giant 
EELR encompasses the bridge. Although we cannot discard 
the possibility that the entire EELR is photoionized by the 
AGN, the spatial variation of the line ratios can be explained 
if the gas in the inner regions (r <24 kpc) is predominantly 
ionized by the AGN while stellar photoionization becomes 
dominant as we move outwards. The star formation might 
have been triggered during the interaction. Whatever ion- 
izes the gas, we find that the properties (sizes, line widths, 
line ratios and luminosities) of the individual knots isolated 
within the LSBN are consistent with nearby HII galaxies. 

Based on the environmental properties (the tidal bridge 
and the companion galaxy "Gl"), we propose that the ori- 
gin of the EELR in this object is tidal debris from a galac- 
tic encounter. The companion must contain a substantial 
amount of metal poor gas, to explain the subsolar EELR 
abundances. The nuclear activity might have been triggered 
in such encounter, rather than in the final stages of a major 
galaxy merger (as the nuclei coalesce). Probably the closest 
approach in the encounter occurred >100 Myr ago, assuming 
a typical 300 km s~^ separation speed, and a true distance 
between the two galaxies of 100 kpc. Such a scenario would 
be consistent with recent imaging results for powerful radio 
galaxies (Ramos Almeida et et. 2010, in prep.). 
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